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Rare-earth ions in dielectric crystals are interesting candidates for storing quantum states of 
photons. A limiting factor on the optical density and thus the conversion efficiency is the distortion 
introduced in the crystal by doping elements of one type into a crystal matrix of another type. 
Here, we investigate the system Pr 3+ :La2(WC>4)3, where the similarity of the ionic radii of Pr and La 
minimizes distortions due to doping. We characterize the praseodymium hyperfine interaction of the 
ground state ( H4) and one excited state ( D2) and determine the spin Hamiltonian parameters by 
numerical analysis of Raman-heterodyne spectra, which were collected for a range of static external 
magnetic field strengths and orientations. On the basis of a crystal field analysis, we discuss the 
physical origin of the experimentally determined quadrupole and Zeeman tensor characteristics. 
We show the potential for quantum memory applications by measuring the spin coherence lifetime 
in a magnetic field that is chosen such that additional magnetic fields do not shift the transition 
frequency in first order. Experimental results demonstrate a spin coherence lifetime of 158 ms - 
almost three orders of magnitude longer than in zero field. 

PACS numbers: 42.50.Md, 76. 30. Kg, 76.70.Hb, 76.60.-k, 03.67.Pp 



> 

(N 
(N 

o 



X 



I. INTRODUCTION 

Rare-earth ion-doped crystals (REIC) have recently 
appeared as promising solid state materials for quantum 
information processing. In the field of quantum comput- 
ing, achieved milestones include controlled phase gatesi 
and single qubit arbitrary rotation^. While these exper- 
imental results were performed on single-qubit and two- 
qubit systems, scalable schemes have also been poposed 3 -. 
In the field of quantum memories, devices able to faith- 
fully store and release photonic quantum states have 
been proposed and implemented. Using several differ- 
ent storage-recall protocols^—, high efficiency^, multiple 
photon storage with large bandwidth2»i£ and entangle- 
ment storag o 11 ' 12 were demonstrated in REICs. These 
results rely on the shielding of the 4f electrons of the 
rare-earth ions by closed shells, which reduces dephas- 
ing by the environment, yielding long coherence lifetimes 
(T2) at liquid helium temperatures. For example, opti- 
cal coherence lifetimes of 4.4 ms have been observed in 
Er 3 +:Y 2 Si0 5 M- 

Even longer lifetimes have been reported for rare-earth 
ion hyperfine transitions and accordingly, the qubit in 
REIC based quantum computing and memories is gen- 
erally defined by selecting two ground state hyperfine 
levels. Optical transitions are used to selectively ad- 
dress qubits or to transfer coherences from the optical 
to the radio-frequency (RF) domain and vice versa. Co- 



herence lifetimes can be extended to 30 s for a ground 
state hyperfine transition of Pr 3+ :Y2SiC>5 at liquid he- 
lium temperature^. This was achieved in two steps: 
first, an external magnetic field was applied to the sample 
in order to decouple one hyperfine transition from mag- 
netic field fluctuations due to host spin flips. As these are 
the main source of dephasing for Pr 3+ hyperfine transi- 
tions, the zero field coherence lifetime of 500 /zs was ex- 
tended in this way to 82 ms^ and later to 860 msii. The 
decoupling was achieved by minimizing the transition en- 
ergy dependence with respect to the magnetic field. This 
condition is referred as ZEFOZ (Zero First Order Zeeman 
shift) transitions^ The coherence lifetime was then fur- 
ther increased by RF decoupling pulses, using a modified 
version^ of the Carr-Purcell sequence originating from 
nuclear magnetic resonance. 

Y2Si05 is the most thoroughly studied host in REIC 
quantum information processing. It combines long co- 
herence life times, favored by its low magnetic moment 
density, mainly due to Y nuclear spins, and high oscil- 
lator strengths. A disadvantage of Y-bascd host ma- 
terial is that doping with Pr 3+ or Eu 3+ , leads to rel- 
atively large inhomogeneous linewidths at high doping 
concentrations, which limits the maximal achievable op- 
tical depth. This is an important concern in high effi- 
ciency quantum memories.— To overcome this limitation 
for Pr 3+ , we proposed a La-based crystal, La 2 (W04)3. 
Pr 3+ substitutes La 3+ in this material, both having very 



similar ionic radii (r ia 3+ = 1.18 A, r Pr s+ =1.14 A)^&. 
Compared to Y 2 Si0 5 (r Y 3+ = 1.02 A}i£, doping stress 
is reduced and the inhomogeneous linewidth is 15 times 
smaller in this compound at high Pr 3+ concentrations. 
However, the magnetic moment of lanthanum (2.78 ^b) 
is much higher than that of Y 3+ (-0.14 /ib) and the 
La2(W04)3 magnetic moment density is 7.5 times higher 
than in Y2Si05. It seems that this should be seriously 
detrimental to coherence lifetimes but we measured a hy- 
pcrfine lifetime of 250 fis by Raman-echoes^, which is 
only smaller by a factor of w 2 compared to the value 
in Y2Si05. This allowed us to measure narrow and ef- 
ficient electromagnetically induced transprency in this 
material^ This result also suggested that REICs that 
could be useful for quantum information processing are 
not limited to the few crystals with very low magnetic 
moment density. However, since applications require T% 
values in the ms range, the techniques described above 
for increasing the coherence lifetime by several orders 
of magnitude should be used. In this paper, we show 
that using a ZEFOZ transition, hypcrfine T 2 can reach 
158 ± 7 ms, corresponding to a 630-fold increase. It is 
therefore possible to strongly reduce the influence of host 
spin flips even in the case of high magnetic moment den- 
sity. 

ZEFOZ transitions appear at specific magnetic fields 
(magnitude and direction), which can only be predicted 
if the system Hamiltonian and all of its parameters are 
known with high precision. In the present system, the 
/ = 5/2 nuclear spin of 141 Pr (100 % abundance) and the 
C\ site symmetry result in a complicated hyperfine struc- 
ture. We therefore used the approach of Ref. [2l|, which 
consists in determining the spin Hamiltonian parameters 
by coherent Raman-scattering before numerically identi- 
fying ZEFOZ transitions. Finally, the coherence lifetimes 
of the hyperfine transitions were measured by optically 
detected Raman-echoes. 



II. MODEL FOR THE HYPERFINE 
INTERACTION 



approximation, the four last terms of Eq. ([T]) can be well 
approximated by a nuclear spin Hamiltonian^ 3 -: 



H n = -g 2 jV 2 B B ■ A • B 

-B ■ {2Ajgjfi B A + gifi N ~E) ■ I 
+1- (P - A 2 7 A) ■ I. 



(2) 



Here, gj is the Lande g-value, [Ib the Bohr magneton 
and B the external magnetic field vector. The A tensor 
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is calculated by second order perturbation theory, a and 
/? denote the coordinate axes and the sum runs over the 
2 J + 1 states of the relevant crystal field multiplet. The 
crystal field levels of interest in the following are the low- 
est of the 3 H4 and *D2 multiplets and are written as 
|0). Other crystal fields levels of each of these multiplets 
are denoted by \n) and E n represents the energy of that 
state. Aj is the hyperfine constant of the J multiplet, gi 
is the nuclear g-factor, /i/v the nuclear magneton, E the 
3x3 unity matrix and / = {I x , I y , I z ) the vector of nuclear 
spin operators. P represents the pure nuclear quadrupole 
tensor. The operator {—A 2 jl ■ A • I) has the same form 
as i" • P • J and is therefore conventionally referred to as 
the pseudoquadrupole interactional. We combine these 
two operators into an effective quadrupole Hamiltonian 
/ • Q • /. In systems with low site symmetry, as in the 
present case, the principal axis systems (PAS) of P and 
A do not coincide in general. The combined operator Q 
therefore has again a different PAS. 

Considering the first term of Eq. ©, we note that it 
does not depend on the nuclear spin. It shifts all nu- 
clear spin levels by the same amount and we therefore 
neglect it in the following. Furthermore, we introduce 
the abbreviation M = (2Ajg,//.isA + gj^ipfE) for the 
effective Zeeman coupling tensor. We thus arrive at the 
following form of the nuclear spin Hamiltonian: 



A good approximation for the Hamiltonian of many 
rare-earth doped compounds is22 



Ho = [Hfi + Ucf] + [Hhf + H q +H 2 



%z 



(1) 



The first two terms, the free ion (including spin-orbit 
coupling) and the crystal field Hamiltonians determine 
the energies of the electronic degrees of freedom. The 
terms in the second bracket, consisting of the hyperfine 
coupling, the nuclear quadrupole coupling, the electronic 
and the nuclear Zeeman Hamiltonian, lift the degeneracy 
of the nuclear spin states. 

The site symmetry (Ci) of our system is low enough 
that the electronic states are nondegenerate. As a result 
of this "quenching" of the electronic angular momentum, 
the electronic Zeeman %z and hyperfine interaction %hf 
contribute only as second order perturbations. In this 



H = B-M-I + I-Q-I. (4) 

The tensors M and Q can be parametrized as^ 
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where the Ri = R(on, /3j,7j) represent rotation ma- 
trices and Euler angles, specifying the orientation be- 
tween M and Q tensors PAS ((x',y', z') and (x",y", z") 
respectively) relative to the laboratory-based reference 
axis system (x, y, z) (see section IIIII) . In general, 
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FIG. 1. Optical transition and hyperfine level structure of 
Pr 3+ :La2(WC>4)3. The order of the energy levels follows from 
previous worki^, whereas the transition frequencies could be 
measured with higher precision utilizing zero field Raman- 
heterodyne scattering within the present work. The arrows 
on the right indicate the range of hyperfine transitions excited 
by RF in the separate experiments. 



the M and Q principal axes are not aligned and ac- 
cordingly the Rq and Rm matrices are not identi- 
cal. The nuclear spin Hamiltonian of a given crys- 
tal field level therefore depends on 11 parameters: 
gx,g v ,gz,aM,f3M,lM,D,E,aQ,f3 Q ,-/Q. 

In zero magnetic field, the quadrupole interaction re- 
sults in a partial lifting of the nuclear spin states degen- 
eracy. The corresponding structures for the 3 H4 and 1 D 2 
levels was determined by holeburning experiment s 19 ! 25 
and are shown in Fig. [I] We label the levels by their pro- 
jections onto the z" principal axes of the Q tensor, noting 
that these states are not eigenstates of the nuclear spin 
Hamiltonian. Since the PAS of the Q tensors of differ- 
ent crystal field levels do not coincide, their quantization 
axes are also different. For the small magnetic fields used 
to determine spin Hamiltonian parameters, the hyperfine 
structure remains close to the zero field one, which allows 
us to identify resonance lines with transitions between 
zero-field states. 



III. EXPERIMENT 

We used Raman-heterodyne scattering (RHS) to ob- 
tain hyperfine spectra of both, the electronic ground state 
and the electronically excited state. RHS is a magneto- 
optic double resonance technique, that requires a laser 
and an RF fiel d 26 ' 27 . The laser light has multiple func- 
tions in the scheme. The first is to prepare a population 
difference between the hyperfine levels by optical pump- 
ing or transfer to auxiliary states. If the RF field is reso- 



nant with the hyperfine transitions, it creates coherences 
between those states, which are transferred to coherences 
in optical transitions by the light. These optical coher- 
ences represent electronic dipoles, which act as a source 
of a new optical field that is shifted by the RF frequency 
with respect to the incident laser frequency - the Raman- 
field. As this field is emitted in the same optical mode as 
the incident light, it can be measured by optical hetero- 
dyne detection, using the transmitted laser beam as the 
local oscillator. 

We used a sample of high optical quality, grown by 
the Czochralski method, containing 0.2% at. Pr 3+ . The 
5x5x5 mm crystal was mounted in an optical cryostat and 
cooled to liquid helium temperatures. La2(WC>4)3 forms 
a monoclinic crystal with a C2/c space group, identical to 
that of Y 2 SiC>5. The La 3+ ions occupy only one crystallo- 
graphic site of C\ symmetry. In each unit cell (containing 
4 formula units), this site appears at 8 positions which 
are related by inversion, translation and Ci symmetries. 
The C2 axes are identical to the C-2 crystal symmetry 
axis, also denoted by b in the following. The C2 sym- 
metry divides the La positions into two groups of 4 ions, 
which behave differently, unless the magnetic field is per- 
pendicular or parallel to the b axis. These two groups 
are called sub-sites in the following. The crystal surfaces 
where polished perpendicular to the (A, Y, Z) principal 
axes of the optical indicatrix. Optical back reflection at 
the Z-surface and mechanical alignment of the X and Y 
surfaces was used to align the crystal along our reference 
(x,y,z) axes, defined by the static magnetic field coils 
(see below). Apart from alignment errors the (A, Y, Z) 
axes should be a replica of (x, y, z), the laser propagating 
along Z and the b axis expected to be closely aligned to 
the laboratory frame y axis. 

A Coherent 899-21 dye laser, further stabilized by 
home-built electronics with respect to intensity and fre- 
quency (linewidth < 20 kHz), served as light source. It 
was tuned to a wavelength of 602.741 nm (vac), corre- 
sponding to the center of the transition involving the 
lowest energy crystal field levels of the 3 H4 and ^2 
multiplets (see Fig. [IJ . Typical powers for the scatter- 
ing/heterodyne light were 0.3 to 2 mW focused, from 
a collimatcd beam of 1.5 mm diameter, with a 300 mm 
lens into the sample. We generated the optical pulses and 
frequency chirps by double pass acousto-optic modulator 
setups. 

The RF fields were applied to the sample by a 10 turn 
6 mm diameter coil. For continuous wave experiments 
(ground state), one side of the coil was terminated by 
a 50 Q load, the other was attached to an RF driver 
running at a power level of 1 W. 

For the exited state spectra, we used a pulsed RHS 
scheme. Here the coil was part of an appropriate tuned 
tank-circuit and typical RF powers were 250 W, resulting 
in maximum signal amplitudes for RF-pulse durations of 
3-4 [is. Figure [3] shows the sequences used for the two 
types of experiments. 

The frequencies for the RF excitation and the shifting 
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The static magnetic field was created by a set of three 
orthogonal Helmholtz coil pairs. They are mounted out- 
side the cryostat and their coil-diameters range from 20 
to 40 cm, providing a homogeneous field over the sample 
volume in their center. With currents of about 10 A, each 
coil pair generates a static magnetic field of about 8 mT. 
To control the field vector a computer control was set up 
for the current sources of the Helmholtz coils. To com- 
pensate non-linearities and drifts, we used a set of three 
orthogonal Hall probes as sensors for a computer-based 
feedback loop. The absolute error of the field compo- 
nents is < 0.06 mT and the relative linear error for the 
static magnetic field is < 0.3%. To minimize the effect 
of small background fields (e.g. earths magnetic field) 
a small compensation field was used, which minimized 
the observed zero-field RHS line splitting and also led 
to almost perfect destructive interferenc e 28 ! 29 . We used 
this compensation field as our zero-field reference in all 
measurements. 

For an optimal determination of the Hamiltonian pa- 
rameters, it is important to sample different strengths 
and orientations of the magnetic field. In our experi- 
ments, we used a spiral on the surface of an ellipsoid^: 



FIG. 2. Pulse sequences. Black lines indicate frequencies and 
grey areas the applied optical/RF powers. 

a) CW (continuous wave) sequence for the ground state 
measurements. A low power laser probe (P p — 0.3 mW) 
and a scanning frequency RF (Prf ~ 1 W, Vi — 7.4 and 
v% = 22.4 MHz for |±±) <-> |±|) or respectively 17.1 and 
32.1 MHz for |±|) o |±|)) were applied to the sample 
to detect the spectra on the frequency encoded time-scale 
(t s — 50 ms). The optimum temperature of the cold finger 
for this scheme was 4.5 K, since still lower temperatures gave 
such slow hyperfine level relaxation rates that it would be 
necessary to repump the hyperfine level population. 

b) Pulsed RHS-sequence for the excited state. Probe and 
erase beam were overlapped in the sample at angle of 0.6°. 
To allow for higher repetition rate the chirped erase laser 
(A = 64 MHz, r e = 10 ms, P e = 8-30 mW due to the fre- 
quency dependence of the AOM) redistributed the popula- 
tions. Initial hyperfine population was created by the probe 
beam (P p = 1.2 mW, r p = 100 fjs) and converted to coher- 
ences by an RF pulse (v RF = 4.94/7.23 MHz (|±|) -H- |±|) 
resp. |±|) o |±| )), Prf = 214/287 W, r RF = 4 /js). For 
optical heterodyne detection the same probe beam was left 
active for an additional time Tdec- The temperature of the 
cold finger for this scheme was 2.4 K. 



of the laser frequency were generated by 48 bit / 300 MHz 
direct digital synthesizers. We controlled the timing of 
the pulses and frequency-chirps by a wordgencrator with 
a resolution of 4 ns. Detection of the heterodyne beat 
signal was accomplished by a 100 MHz balanced photo 
receiver (Femto HCA-S), a phase sensitive quadrature- 
detection demodulation scheme, appropriate analog and 
digital filters and a digital oscilloscope. 



B(t) 



B x Vl 

ByVl 



- t 2 ■ cos (67rf) 

- t 2 ■ sin (6irt) 
B z -t 



(7) 



Here, we use 



t=-l + (JV-l) 



Ntnt ~ 1 



N = 1,2,...,N U 



to represent the discrete coordinate along the tra- 
jectory. For the ground state series, we measured 
Ntot = 101 orientations, with magnetic field amplitudes 
[B x ,B y ,B z ] = [7,9,8] mT and for the excited state we 
used Ntot = 251 and B x = B y = B z = 6.5 mT. Figure [3] 
shows some typical experimental spectra for the ground- 
and excited state. 



IV. THE FITTING PROCEDURE 

In the case of Pr 3+ :La2(WC>4)3, 11 spin Hamiltonian 
parameters are necessary (D, E, gx,y,z, (a, /3,j)q/m) to 
compute the line positions of excited or ground state 
spectra for an arbitrary external magnetic field. These 
calculations have to take into account the magnetically 
non-equivalent Pr 3+ sites, which arc related by the crys- 
tal Gi/b axis, corresponding to the Indirection of the 
crystal's optical indicatrix axis system (X, Y, Z). As this 
system may not be perfectly aligned to the laboratory 
reference system (x, y, z) (see section UTI]) . (a,j3)c 2 an- 
gles give the Ci orientation with respect to the x and z 
axes. This gives a total of 13 parameters that are to be 
extracted from the recorded series of spectra. 
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FIG. 3. (Color online) Representative ground state CW RHS and exited state pulsed RHS spectra. Both spectra from the 
ground state (\g,i «-> j)) are recorded at N = 15, corresponding to B = (2.60, —5.27, —5.76) mT. For the excited state spectra, 
we used N — 211, corresponding to B = (4.62, 1.19,4.42) mT. The normalization is relative to the largest line from ground 
or excited state spectra, respectively. The ground state spectra shown here resolve all 8+8 possible RHS transitions, while 
in the excited state spectra not all lines are resolved. The histograms on the right show the distributions of fitted full width 
at half maximum (FWHM) RHS-linewidths for all recorded data (see Section [IV A[) . plotted separately for the two hyperfine 
transitions (\g/e,i o j), i/j = ±5/ ± | or ±|/ ± §). For both, the ground and the excited state 
are bigger (see Section IVCj) . The mean linewidths are ki| <-> ±7 

o ±|) w 18.3 kHz, je,±| o ±f ) « 26.3 kHz, |e, all) « 22.3 kHz 
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301 kHz, |g, all) « 196 kHz 
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B. Fitting of Hamiltonian parameters 



While it is possible to calculate directly the spec- 
tra from the Hamiltonian, the resulting amplitudes and 
phases do not agree with the experimentally observed 
values, since they depend strongly on experimental de- 
tails, such as cable lengths, tuning curves etc. By far the 
most precise parameters of the spectra are the positions 
of the resonance lines. We therefore extracted these in 
a first step and used only the resonance frequencies as 
input to the actual fitting procedure. 

For the extraction of the line positions, we fitted the 
absolute- value representation of the measured spectra to 
a series of gaussian-shaped resonance lines. To discrimi- 
nate between real lines and noise or artifacts, we consid- 
ered additional parameters like peak height, line width 
and noise level. For most of the recorded data the pro- 
cedure worked satisfyingly, even for partially overlapping 
resonance lines. The positions of lines that had very small 
amplitudes or could not be fitted well due to strong over- 
lapping were identified visually or omitted. 



After determining the resonance line positions, we fit- 
ted the Hamiltonian parameters such that the resulting 
line positions agreed with the experimental data. Due 
to the large number of parameters and their complicated 
intcrdependency, gradient-based algorithms are not effi- 
cient, as they tend to stick to local minima. In our case, 
the combination of first running a probabilistic and then 
a direct search yielded the best convergence to the global 
minimum. As in the work of Longdel l 16 i 21 we used simu- 
lated annealing^ for the first step. When the preliminary 
results of this step appeared to converge towards a global 
minimum, we switched to a pattern search algorithm^!, 
which polls meshes of adjacent points to find a better 
minimum. This reduced the necessity for an excessively 
long search at low temperatures in the simulated anneal- 
ing procedure and represents an additional check for the 
quality of the minimum. 



The implemented algorithm minimizes the root mean 
square (RMS) deviation between the measured and the 



calculated line positions of all spectra, 
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For each of the N tot spectra Qn line pairs, consisting 
of the measured frequencies v^ p and the corresponding 
calculated frequencies v^ N were identified. The index i 
runs over all lines in a single spectrum. Accordingly all 
possible eigenvalue differences {v^ N ) for the given Haniil- 
tonian, derived from Eq. (J4J using the given set of pa- 
rameters of the k th iteration and the N th magnetic field 
orientation (Eq. ©), were calculated. Not all possible 
lines are always resolved or visible in the spectra (see 
e.g. Fig. [3]). Therefore, we matched theoretical to ex- 
perimental lines by minimizing the absolute deviations 
for this spectrum; each experimental line was associated 
only with a single theoretical line. Lines that could not 
be assigned uniquely were not considered in that step of 
the fitting procedure. This mapping could change during 
the fitting procedure, but in the final steps of the search, 
the procedure resulted in an unambiguous assignment. 

We initialized the simulated annealing algorithm with 
a random guess for the parameters. At each following 
iteration one of the Hamiltonian parameters was ran- 
domly chosen and varied in an interval determined by 
the current temperature. If the variation resulted in a 
new minimum value of the RMS error (/^^f ), the new 
point was always accepted. If the variation was higher 
than the previous minimum, the algorithm calculated the 
Boltzman factor 



p k = exp 



I fRMS _ fR,MS\ 
\Jk Jbest ) 

kn, ■ Tu, 



(8) 



and compared it to a random number p r e [0, 1]. The new 
point was accepted for pk > p r and rejected otherwise. 
The RMS misfits f^ MS and f^ e f t s represent the devia- 
tion of the current iteration and the best observed devia- 
tion of all previous iterations, respectively. The temper- 
atures Tk, are defined in the individual l-th parameter 
units. Conversion factors fcs, ensured correct units and 
parameter independent scaling of the denominator in Eq. 
(|8]). All temperatures were lowered continuously as the 
fitting improved. This was done for all Tj., in the same 
way, so that for simplicity we may speak in terms of a 
global temperature and omit the Z-index where it is not 
necessary. At high temperatures the algorithm sampled 
a large parameter range. Enabling it to resettle in states 
of slightly higher misfits enables the algorithm to trace a 
large search-space and to avoid getting stuck in a local 
minima at the same time. As the temperature was gradu- 
ally lowered, the parameter values became more confined 
in the vicinity of the global minimum. 

The whole procedure was implemented as a MATLAB 
program, utilizing the Global Optimization Toolbox 
functions patternsearch and simulannealbnd. 



We chose the initial temperature To such that the cor- 
responding changes of the line positions were « 10 times 
the inhomogencous line widths. Typically k max « 10 6 
iterations for the simulated annealing lead to reliable re- 
sults. Within the first 80% of the iterations the temper- 
ature was gradually lowered as T& « To (1 — k/k max ) . 
Then Tk was left constant at a value corresponding to 
a frequency uncertainty of a fraction of the typical line- 
widths (a few kHz). 

During the whole fit, the parameters were constrained. 
As values for D and E can be derived from hole- 
burning experiment o 19 ' 25 their boundaries were restricted 
to ±10% of the literature values. For the gyromagnetic 
ratios, the boundaries were chosen to ±100% of the ex- 
pected values of -(10-100) MHz/T— observed in similar 



systems r 



The Euler angles were allowed to vary over 



the whole definition range. The probability for choosing 
a specific parameter for the variation was proportional to 
the (relative) size of its boundaries. As indicated above, 
the new value of this parameter was chosen by adding a 
random value out of the interval ±Tfc, . If the new value 
did not lie within the boundaries, the process was re- 
peated. The final direct search was set up with the same 
boundaries and typically terminated after 5 • 10 itera- 
tions. 

The orthogonal axes of the B-field coils define the refer- 
ence coordinate system (x,y, z). Euler-angles and trans- 
formations refer to this basis and are given in "zyz"- 
convention^ (see. Eq. (|A.1[) for details). As indicated 
before, we did not constrain the relative orientation of 
the quadrupole Q- and Zeeman M-tensor. To avoid am- 
biguous results, we fitted all measured spectra of a sin- 
gle electronic state simultaneously. Due to the two non- 
equivalent sites this results in a maximum of 16 lines per 
spectrum. 

Since the assignment of the resonance lines to the two 
sites is not known, we had to fit both sites simultaneously. 
Instead of fitting 2 ■ 11 — 22 parameters, we used the fact 
that they are related by a C<i rotation. Using Eq. ([4]) for 
site 1, we write the Hamiltonian for site 2 as 

n 2 = b ■ {RcMiRl 2 ) ■ i + 1- (Rc 2 QiRc 2 ) ■ I 

with 

Rc 2 — Re ' Rk ' Re } 

Rc = R(a C2 ,Pc 2 ,0), R W = R(180°,0,0). 

The angles ac 2 and Pc 2 correspond to the spherical co- 
ordinates of the Ci axis in the laboratory system. 

V. RESULTS 

The underlying symmetry of the crystal field and the 
structure of the spin Hamiltonian cause some ambiguity 



TABLE I. Ground state spin Hamiltonian parameters 



(SNRi) for all contributing RHS-lines: 



parameter 


value 


fit error 


unit 


D 


-6.3114 


0.0027 


MHz 


E 


-0.8915 


0.0021 


MHz 


Q.Q 


20.4 


3.3 


deg. 


PQ 


147.7 


1.4 


deg. 


1Q 


10.2 


1.4 


deg. 


9x 


-51.7 


3.6 


MHz/T 


<hi 


-23.5 


1.1 


MHz/T 


9* 


-146.97 


0.75 


MHz/T 


Q-M 


30.1 


3.8 


deg. 


Pm 


146.59 


0.55 


deg. 


JM 


13.09 


0.69 


deg. 


0iO 2 


88.34 


0.47 


deg. 


Pc 2 


92.45 


0.31 


deg. 



if only RHS spectra are used to determine the Hamil- 
tonian parameters. The thesis of J. Longdell 3 ^ provides 
a detailed review of relevant symmetries that make it 
impossible to unambiguously determine all Hamiltonian 
parameters from RHS spectra alone. Important for our 
investigation is the fact that the RHS spectra do not 
depend on the signs of D, E and the gyromagnctic fac- 
tors g x , g y and g z . In addition, different sets of Eulcr 
angles correspond to the same tensor orientations. As 
a consequence, different runs with random initial values 
lead to apparently different solutions. We checked that 
these solutions arc related by the symmetry operations 
mentioned above and verified thus that we really found 
a unique global minimum. 



i L 



For the ground state data we found v a = 1.4 kHz. Ac- 
cording to this we chose the fixed temperature T ai so 
that a single parameter change from its optimum value 
by T ai resulted in an increase of the RMS deviation by 
v a . After 2 • 10 6 iterations the histograms of the accepted 
parameters all showed a gaussian shape, whose lcr-width 
are given as fit error in Table |TJ 

Apart from the statistical error, we also consider sys- 
tematic errors. The most important contribution is due 
to the calibration error of the magnetic field. We estimate 
its precision to ~ 0.65%, which translates to the same 
fractional uncertainty of the gyromagnctic ratios g x , g y 
and g z . As the parameters are given in the laboratory- 
fixed reference frame (x, y, z) a misalignment of the crys- 
tal does not contribute to the error but is expressed by 
the ac 2 and /3c 2 values. They determine the orienta- 
tion of the Ci axis and therefore also that of the op- 
tical indicatrix (X, Y, Z) in our laboratory frame. The 
only systematic contribution in the angles arise from non- 
orthogonality of the coils, which is < 1°. The uncertainty 
of the alignment of the crystal relative to our reference 
(x, y, z) and that of the crystal surfaces to the optical in- 
dicatrix (X, Y, Z) results in an error of « 5° for the angles 
seen relative to the crystal axis system. The frequency 
scan was generated by direct digital synthesis, resulting 
in negligible uncertainty in the frequency scale. 



B. Excited state 



A. Electronic ground state 

Figure [1] shows the experimental data for N to t = 101 
different external magnetic fields. Several fit trails reli- 
ably led to the parameters shown in Table U and repre- 
sented by the solid lines in Fig. 01 

With these parameters, the RMS deviation between all 
accounted line positions and the fit is f=s 32 kHz, signifi- 
cantly smaller than the average linewidth of the ground 
state RHS lines of « 196 kHz (see Fig. [3]), indicating 
that it is dominated by statistical error. At the end of 
the fitting procedure, L =1218 of the total 1221 experi- 
mental lines from N to t = 101 spectra could be assigned 
to calculated resonance line positions. To estimate the 
uncertainty of the fitted parameters, we sampled the pa- 
rameter space in the vicinity of the global minimum by 
repeating the probabilistic part of the fitting procedure 
again using a fixed, low temperature. Such a procedure 
can be shown to be rigorous if the only source of error is 
gaussian noise in the line positions! 21 ' 34 As an estimation 
for the noise in the line positions we used the mean ratio 
of fitted line widths Oi to the individual signal to noise 



For the excited state we used the same procedure as for 
the ground state. With the optimal fit parameters, we 
found an RMS deviation of 3.1 kHz between theoretical 
and experimental frequency values, using L = 2345 of the 
2353 measured lines in N to t = 251 spectra. Compared to 
the mean experimental FWHM of 22.3 kHz, the RMS 
deviation is even better than for the ground state. We 
mainly attribute this to the higher quality of pulsed RHS 
spectra, with fewer line shape artifacts. Figure [5] and 
Table ITT1 show the results. 

For the determination of the fit errors, we used the 
same procedure as for the ground state. The average 
uncertainty of the line positions was v a = 136 Hz. The 
systematic errors are again dominated by the calibration 
error of the magnetic field. Although the crystal was 
remounted between the excited and ground state exper- 
iments, the resulting orientation of the C^ axis agrees 
between the two data sets within less than one degree, 
which is less than our alignment accuracy. As the Zee- 
man tensor is almost axially symmetric, g x ~ g Zl its 
orientation relative to the quadrupole tensor or to the 
ground state tensor orientations can not be completely 
determined. 
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FIG. 4. Ground state |±§) -H- |±|) and |±|) o |±|) liyperfine transition CW RHS spectra. The solid lines represent the 
fit results and the shaded background the absolute value of the experimental spectra. For each field orientation, the spectrum 
was normalized to the maximum signal amplitude. 
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5. Excited state l±|) 



|±|) and |±|) <-> |±|) hyperfine transition pulsed RHS spectra. The solid lines represent the 



fit results and the shaded background the absolute value of the experimental spectra. For each field orientation, the spectrum 
was normalized to the maximum signal amplitude. 



TABLE II. Excited state spin Hamiltonian parameters 



parameter 



value 



fit error 



unit 



D 


1.90705 


0.00023 


MHz 


E 


0.35665 


0.00014 


MHz 


Q-Q 


-18.51 


0.71 


deg. 


PQ 


73.83 


0.48 


deg. 


1Q 


-84.22 


0.37 


deg. 


9x 


-17.22 


0.27 


MHz/T 


flu 


-14.39 


0.10 


MHz/T 


(J-- 


-18.37 


0.14 


MHz/T 


OtM 


-23.7 


2.4 


deg. 


Pm 


88.5 


4.1 


deg. 


JM 


-80.1 


1.8 


deg. 


OtC-, 


88.63 


0.25 


deg. 


Pc 2 


92.69 


0.24 


deg. 



TABLE III. Crystal field parameters £>,,. calculated Q and M 
tensor principal values and second order hyperfine interaction 
parameters. 



B i:j 


(cm" 1 ) 




ground state 


excited state 


unit 


B20 


375 


D 


-6.1 


2.0 


MHz 


B22 


-93 


E 


-0.63 


0.30 


MHz 


B40 


768 


9xa 


-32 


-18 


MHz/T 


B 42 


445 


9Va 


-22 


-4 


MHz/T 


B44 


1027 


9zc 


-151 


-18 


MHz/T 


Beo 


267 


Aj 


937 


697 


MHz 


B@2 


-402 


9J 


0.81 


1.03 




B@4 


-61 


l±-x c x c 


0.0280 


0.0066 


cm 


B ee 


-52 


^-VcVc 


0.0140 


-0.0090 


cm 






A 2c2c 


0.2000 


0.0069 


cm 



C. Discussion 

Examining Tables U and [Til it appears that the prin- 
cipal values for the Q and M tensors are similar to 
those found in other hosts like Y2Si0^i, LaF3 or 
YAIO^. Especially, the ground state gyromagnetic ten- 
sor is anisotropic with one large component, in contrast 
to the excited state which also exhibits smaller values. 
To get some insight into these properties, we compared 
these results with calculations derived from crystal field 
calculations. 

In a previous work 3 ^, we found that quadrupolar D 
and E values for ground and excited states (see Eq. (|6])) 
could be very well reproduced starting from electronic 
wavefunctions obtained by a crystal field analysis. The 
latter was done assuming a C2V site symmetry, which is 
higher than the actual one (Ci). In this higher symmetry, 
the Q tensors for different crystal field levels are colinear, 
in clear contradiction with our results. Nevertheless, it 
seems that the additional crystal field parameters of C\ 
symmetry have little effects on the Q principal values. 
In the following, we present the calculations of the M 
tensor principal values. 

In C2V orthorhombic symmetry, the spin Hamiltonian 
of Eq. (J4]), expressed in the (x c , y c , z c ) crystal field axes, 
reads: 



H" = J2 B ^ h 



-D(ll 



^Wi-4 



where the gi are related to the A tensor by (see Eq. ©) 



g t = -2AjgjHBKi - gi^N- 



(10) 



The A tensor is given by Eq. (O and can be calculated 
from the electronic wavefunctions. The latter were found 
using a free ion and crystal field Hamiltonian whose pa- 
rameters were fitted to experimentally determined crys- 
tal field levels 3 -^. In this calculation, we use arbitrary 



permutations of the (x c ,y c ,z c ) axes. This results in dif- 
ferent sets of E,D values which give the same hyperfine 
energy levels. We subsequently fix the choice of the axis 
system such that the convention < 3E/D = rj < 1 
|36| is fulfilled, thereby resolving ambiguous sets of pa- 
rameters, such as (in MHz) D = 2.082, £ = -3.3445 
and D = -6.0578, E = -0.6313, which describe identi- 
cal ground state zero field hyperfine structures and corre- 
spond to an exchange of x c with z c . With this convention, 
we fix the permutation of the axes and thus the values for 
D and E for both electronic states. The crystal field pa- 
rameters we used and the corresponding E and D values 
for the ground and excited states are listed in Table ITTTF ' 
The electronic wavefunction of the excited 1 D2(|1)) level 
is used instead of that of 1 D2(|0)) to take into account 
a wrong ordering in the calculated crystal field levels of 
this multiplet^^. 

Comparing Tables HI and Hll with Table ITTTl shows that a 
reasonable agreement is found between experimental and 
calculated principal values of ground and excited state 
M tensors. Again this suggests that additional parame- 
ters appearing in calculations using C\ symmetry mainly 
determine the relative orientation between the different 
tensors. Calculated values especially reproduce two fea- 
tures mentioned above: the very large value of g z for the 
ground state and the smaller gi values for the excited 
state compared to the ground state. A qualitative un- 
derstanding of these properties can be obtained from the 
crystal field analysis by looking at the different factors en- 
tering in Eq. ([TO]) . They are summarized in Table Unl We 
first note that the isotropic and crystal field independent 
nuclear Zeeman contribution to gi equals -12.2 MHz/T. 
Differences in gi values are mainly linked to the pscudo- 
quadrupolc A tensors since the products Ajgj vary by 
only 4% between ground and excited states. The pseu- 
doquadrupolc tensors involve the J matrix elements and 
the energy differences appearing as denominators in Eq. 
(J3]). We first discuss the ground state case. The elec- 
tronic wavefunction of interest (the lowest energy crystal 
field level) has the following form: 

|0) = -0.62 \ 3 H±, -4) - 0.62 | 3 ff 4 , 4) - 0.4 | 3 ff 4 , 0) 
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where brackets on the right hand side are written as 
\ 2S+1 Lj, Mj) and only terms with a coefficient larger 
than 0.15 have been kept. The larger J Zc matrix element 



is found between 
only composed of 



0) and |1), since the latter is nearly 
3 iJ 4 ,±4) states. The |(0|J z Jl)| ma- 
trix element equals 3.6 close to the maximum value of 
|( 3 -ff4,±4| J 2c | 3 i?4,±4)| = 4. Moreover, this large ma- 
trix element is found for levels close in energy (65 cm -1 
Ref. l35l ). resulting in a large A ZcZc . On the other hand, 
|0) couples to crystal field levels containing | 3 7l4,±3), 
| 3 iJ4,±l) by J Xc or J Vc operators. The corresponding 
matrix elements do not exceed 2.4 in absolute value. 
As expected, this is close to the average value of ma- 
trix elements of the form ( 3 i?4, ±4| Ji | 3 i?4,±3) and 
( 3 -ff4,±l| Ji | 3 i?4,±2) (where i = x c or y c ), which is at 
most 1.9. The levels with the largest matrix elements are 
located at high energies (#3 = 143 and E 5 = 349 cm -1 
for Jy c and J Xc respectively), resulting in low A XcXc 



A 



and 
and 



y c y c . This in turn explains the small values of g x 
g yc compared to g Zo . 

A similar analysis can be performed for the X D 2 excited 
state. The level of interest is |1) because the crystal field 
calculation inverts levels |0) and |1) as mentioned above. 
The latter is found to be equal to: 

|1> =0.67\ 1 D 27 -2)-0.67\ 1 D 2} 2) 

with the same convention as above. This state gives a 
J Zc matrix element equals to 2 with the state |4), lo- 
cated 441 cm -1 higher than |1). Maximum average val- 
ues for matrix elements of J Xc and J y< , can be estimated as 
above for levels containing \ 1 D 2 ± l) states, resulting in 
KOlJaJl)! « |(l|J Xo |2)| ps 1. The corresponding ener- 
gies are E — E\ = —82 cm -1 and E 2 — E\ = 113 cm -1 . 
The combination of matrix elements and energy differ- 
ences result in smaller values for An (i — x c , y c , z c ) com- 
pared to the ground state. This can also partly explain 
the isotropy of the excited state gi values, which are closer 
to the nuclear Zeeman contribution. As pointed out 
above, several Pr 3+ doped compounds exhibit the same 
behavior so that the discussion given above could also be 
applied to them. 

We now turn to the principal axes of the spin Hamil- 
tonian tensors. As a further test of the (oq,I3q,jq) 
parameters determined from the RHS experiments, we 
compared the relative oscillator strengths obtained from 
zero field spectral tayloring experiments^ 9 - with calcula- 
tions. The oscillator strengths are assumed to be pro- 
portional to the square of the overlap of the nuclear 
wavcfunctions2£, the latter being given by the ground 
and excited state Hamiltonians. This assumption is rea- 
sonable since the hyperfine interactions are a small per- 
turbation to the electronic wavefunctions. The results 
are gathered in Table [IVJ A good agreement is found, 
showing that indeed the orientation of the quadrupole 
tensors was determined correctly. In Pr 3+ :Y2SiC>5, sig- 
nificant discrepancies were found between calculated and 
experimental value s 39 ! 40 . This was tentatively attributed 



TABLE IV. Experimental and calculated relative optical os- 
cillator strengths between 3 H4 and X D2 hyperfine levels. Rows 
correspond to transitions starting from the ground state hy- 
perfine levels and columns correspond to transitions to differ- 
ent excited state hyperfine levels (see Fig. [TJ . Experimental 
data from Ref. HiA 



\e,±h) 



■±§> 



e ±- 



(s,±i| 

<5-±|| 



exp 
cal 



exp 

cal 



exp 

cal 



0.09 ±0.01 
0.08 ±0.01 

0.33 ±0.01 
0.31 ±0.02 

0.55 ±0.01 
0.60 ±0.02 



0.28 ±0.01 
0.24 ±0.02 

0.39 ±0.01 
0.45 ±0.02 

0.36 ±0.01 
0.31 ±0.02 



0.63 ±0.01 
0.67 ±0.02 

0.28 ±0.02 
0.24 ±0.02 

0.09 ±0.01 
0.09 ±0.01 



to additional selection rules due to superhyperfine cou- 
pling with Y ions. In our case, it seems that although 
superhyperfine coupling may also be observed (see Sec. 
IV Dl) , relative optical transition matrix elements can still 
be determined from the overlap of the nuclear wavefunc- 
tions. 

Hyperfine transition lincwidths were also determined 
during the fit procedure (Fig. 3). The data show 
that the transitions with the larger splittings also show 
the larger linewidths. For example, the ground state 
|±|) O |±|) transitions at 24.44 MHz have an average 
linewidth of 301 kHz (see Fig. 3, caption) whereas the 
I ±|) «-> |±|) transitions at 14.87 MHz have a linewidth 
of only 105 kHz. To explain this, we first consider that 
the used fields arc small enough, so that the observed 
lincwidths are similar to those obtained at zero field. 
Moreover, we approximate the spin Hamiltonian by set- 
ting E = in Eq. © so that H" = D(I* c - J(7+l)/3). In 
this case, the 1 1) <-)• ||) transition energy is \2D\ and the 
||) f-j- ||) is |4D|. Crystal field variations from one ion 
position to an other correspond to a distribution of crys- 
tal field parameters and therefore of the D parameter. 
The hyperfine linewidths in the excited and ground state 
should then be proportional to the transition energies. 
This is qualitatively in agreement with the experimen- 
tal values. The excited state linewidths are also smaller 
than the ground state ones, which suggests that the D 
distribution width is also proportional to D. 



D. Experimental verification of a ZEFOZ transition 

As mentioned earlier the coherence times for ZEFOZ 
transitions arc expected to be much longer than at zero 
or arbitrary magnetic field. Up to date this was demon- 
strated experimentally only for p r 3 + : Y 2 Si05J 14 ' 15 To 
proof the usefulness of the ZEFOZ technique for other 
compounds and also to verify our hyperfine characteri- 
zation we present here experimental data of a ZEFOZ 
transition of Pr 3+ :La2(W04)3 in the following. Using 
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our parameterization of the spin Hamiltonian we sought 
for magnetic field configurations and transitions that sat- 
isfy the ZEFOZ conditions^: 



S'iB. 



opt) 



dvjjBopt) dvi(B opt ) 



9B X 



f)B u 



dvi(B opt ) 
dB z 



0. 



. (11) 

We identified such points by numerical minimization of 
|S r/ (i?)| for all transitions v± within a static magnetic 
field grid. In this way we found several ZEFOZ posi- 
tions where one transition satisfies Eq. (|11|) and further 
showing low curvature, e.g. small second order coeffi- 
cients 



SJi(B) = 



d 2 n{B) 



dBjdBk 



(12) 



An identified (hyperfine ground state) ZEFOZ transition 
at v A = 12.6 MHz and B opt = (57.5, 4.0, -36.1) mT was 
experimentally explored. The setup at TU Dortmund de- 
scribed in Section lnll was not designed for magnetic fields 
of more than 12 mT per axis. Experiments exploring 
the ZEFOZ point we carried out in Lund. This offered 
the opportunity to experimentally verify the Hamilto- 
nian parameters and predicted ZEFOZ points in an in- 
dependent laboratory. The static magnetic field vector 
was provided by a set of three orthogonal superconduct- 
ing coils, the y-coil being part of an Oxford Spectromag 
cryogenic 7 T magnet assembly with 0.1 mT resolution. 
The homebuilt x and z coils could generate fields of a 
few 100 mT and were controlled by 16 bit DAC. To fit 
into the homogeneous region of the coils we had to cut a 
5x5x1 mm piece from the sample that was used for the 
characterization. Due to this and the construction of the 
sample holder, we could only align the optical indicatrix 
with high precision along the z-axis (laser direction) of 
the coil frame. The other axis orientations were only 
known with a precision of about 10 degrees. To find 
the ZEFOZ point experimentally we had to consider this 
misalignment, the accuracy of the Hamiltonian parame- 
ters and the calibration of the coils. Therefore in a first 
step we adjusted the magnetic field to get a good over- 
lap between observed CW RHS spectra and the calcu- 
lated line positions, that follow from Table U and the ZE- 
FOZ field. As the transition of interest shows very small 
frequency changes when being close to the desired field 
vector, we utilized a T 2 measurement by Raman-cchos, 
induced by two RF pulses (P as 3 W, t p ps 25/50 /lis, 
pulses along z-axis), in a second step. Thus we fine tuned 
the magnetic field components for maximum Raman-echo 
signal at long echo times (RF pulse separations). Fig- 
ure [S] shows the longest-lived Raman-echo decay curves 
we could achieve. These demonstrate hyperfine coher- 
ence times T2{B op t) of up to 158 ± 7 ms, representing 
a 630 fold increase compared to the zero magnetic field 
situation™. The decay curves at magnetic fields slightly 
detuned from the ZEFOZ point show slow modulations 
with a frequency of vm = 24.5 ± 1.7 Hz. This could be 
due to a supcrhypcrfine interaction with La nuclei, but 
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FIG. 6. (Color online) Raman-echo decays at ZEFOZ point 
and with magnetic detunings of -0.2 and -0.5 mT for the z- 
component. 



a clear explanation is lacking at the present time. This 
point will be investigated in further experiments. To fit 
decay times we used the function 



f(t) = (A + M cos (2iru M t)) exp 



where A is the exponentially decaying part, vm is the 
modulation frequency, T 2 the decay time and c an off- 
set. As we moved the magnetic field away from the 
ZEFOZ point B opt by —0.2 mT in the z-component 
this resulted in a decrease of the coherence time to 
T 2 (-0.2 mT) = 133 ± 16 ms and a shift by -0.5 mT 
resulted in T 2 (-0.5 mT) = 97 ± 19 ms. 

In the following we will use this measured T 2 (B) val- 
ues to estimate the magnetic field fluctuation AB at the 
Pr 3+ site. These fluctuations cause frequency shifts and 
thereby broaden the hyperfine transition. For this pur- 
pose, we expand the hyperfine transition frequency Vi on 
a deviation -B D // from a given field B 



S =Vi{B) + slB ofs 



.s-2 



B 



off 



If the reference field B fulfills the ZEFOZ condition, the 
first derivative vanishes, si = and the frequency shift 
due to the fluctuations is 



A, = f(Ai?) 2 . 



(13) 



Since the fluctuations only generate positive frequency 
shifts (for s2 > 0), they add a line broadening of Az//2. 
Using the assumption that the line broadening is entirely 
due to magnetic field fluctuations^, we obtain a decay 
rate 



1 2 



»2 



(ABf 
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With the experimentally obtained value of T 2 = 158 ms 
and s2 as 12 kHz/mT 2 , calculated from the maximum 
eigenvalue of the derivative matrix Eq. (fT2|) . using the 
parameters from Table HI we thus estimate the magnetic 
field fluctuations as AB « 46 //T. This value is of the 
same order of magnitude as that found for the ZEFOZ 
points in Pr 3+ :Y 2 Si05, the only ones experimentally in- 
vestigated up to date, where s2 pa 3-6 kHz/mT 2 and 
AB = 14 /i T 14 ' 16 were found. The zero field relaxation 
times in Pr 3+ :La 2 (W04)3 an d Pr 3+ :Y 2 Si05 are also com- 
parable (w 250 /is vs. w 500 ^s). 

We now analyze the dependence of the relaxation times 
on the magnetic field when the deviation is large com- 
pared to the amplitude of the fluctuations, B ff 3> AB. 
This changes Eq. (JTHJ) to 



s2 
T 



--^(B 2 off + 2ABB off) 



The first term describes the line shift, the second a line 
broadening. Using the value AB pa 46 /iT for the mag- 
netic field fluctuations, we expect that the resulting line 
broadening for a field change of B // =0.2 mT is « 9 ms 
and for -B Q // = 0.5 mT ps 3.6 ms. These values are sig- 
nificantly shorter than the experimental values. 

The most likely explanation for this discrepancy is that 
the relaxation at our reference field is not entirely due 
to magnetic field fluctuations and that the reference field 
does not exactly fulfill the ZEFOZ condition. Both effects 
lead to additional contributions to the dephasing rate. 
We therefore write the total dephasing rate as 



t.; 



= T^ + slAB 



s2ABB off . 



(14) 



Here, T 2 _ describes those contributions that are not due 
to magnetic field fluctuations, such as phonons, while 
si is the first derivative of the transition frequency with 
respect to the magnetic field change. Both contributions 
to the dephasing rate are independent of the magnetic 
field offset -Bo// and therefore not distinguishable in the 
available experimental data. 

Using now s2 ps 8.2 kHz/mT 2 , corresponding to the 
projection of Eq. (TT2"]) into the the z-direction, we use 
Eq. (fT4"|) to estimate the magnetic field fluctuations. The 
result of a linear fit of T^ 1 vs. -Bo// yields AB pa 1 /iT. 
Applying the same analysis to the Pr 3+ Y 2 SiOs data (es- 
timated from Figure 2 in Ref. KL4I ) gives AB pa 1.3 /iT. 
This also reduces the estimate for AB by approximately 
one order of magnitude compared to the analysis where 
the minimum dephasing rate is assumed to originate 
entirely from the quadratic term of the magnetic field 
fluctuations^. 

To get a better estimate for AB further measurements 
would be required. The phononic contributions to T 2 
could be determined from measurements at different tem- 
peratures, as si and s2 are independent of this param- 
eter. Quantitative measurements of T 2 at several devia- 



tions B ff in three directions could help to estimate the 
numerical value of the second term of Eq. I|14p. 
VI. CONCLUSION 

We characterized the hyperfine interaction of 
praseodymium doped into La 2 (W04)3 for the elec- 
tronic ground state and one electronically excited state. 
We described in detail the experimental and numer- 
ical methods to reliably derive the spin Hamiltonian 
parameters. The relative oscillator strengths between 
the 3 H4 and 1 D 2 hyperfine levels derived from our data 
are in good agreement with those measured in earlier 
experimental worki^. We indicated physical reasons 
for the experimentally found tensor orientations and 
their principal axis values by a crystal field analysis 
discussion. Further the calculated tensor values fol- 
lowing from this analysis are in reasonable agreement 
with the values derived from the experimental data, 
backing up our results and showing the usefulness of 
such analysis techniques. The full characterization 
enabled us to calculate all transition frequencies for 
arbitrary magnetic fields. Using this, we could predict 
the magnetic field value at which a ZEFOZ transition 
occurs. We verified this condition experimentally in a 
second laboratory. Besides the experimental verification 
of the hyperfine characterization, we determined the 
order of magnitude of the magnetic fluctuations at the 
Pr 3+ -sitc (AB < 10 /zT) and investigated the coherence 
properties of the ZEFOZ transition. The latter showed 
characteristics similar to superhyperfine interaction 
and a up to 630-fold increase of the coherence lifetime 
compared to zero field. The demonstrated spin lifetime 
of 158 ± 7 ms and the relatively low second order 
Zeeman coefficient show that even crystal systems with 
high magnetic moment density can have a high potential 
for quantum memory and information applications. 
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Appendix: Euler angles conventions 

For the tensor transformations we used the common 
"zyz" -convention^ 3 - with right-handed coordinate sys- 
tems. The transformation matrices R, e.g. in Eq. [5] and 
[SI are given by: 
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cos a sin a 
— sin a cos a | , B 
1 



cos ft — sin ft 

1 
sin ft cos ft 



cos 7 sin 7 

— sin 7 cos 7 

1 



i?(a )y 8,7) = C-B-A = 



sin a sin 7 + cos a cos /3 cos 7 cos a sin 7 + cos ft cos 7 sin a — cos 7 sin ft 
cos 7 sin a — cos a cos ft sin 7 cos a cos 7 — cos /3 sin a sin 7 sin ft sin 7 
cos a sin ft sin a sin ft cos /? 



(A.l) 



r 



with a, 7 e [— 7r,7r] and /3 e [0,7r]. 
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